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TEMPERATURE DISTRIBUTION IN INTERNALLY HEATED WALLS OF HEAT 
EXCHANGERS COMPOSED OF NONCIRCULAR FLOW PASSAGES 1 

By E. R. G-. Eckebt and Georges M. Low 


SUMMARY 

In ike walls of heat exchangers composed of noncircular 
passages, the temperature vanes in the circumferential direction 
because of local variations of the heat-transfer coefficients . A 
prediction of the magnitude of this variation is necessary in order 
to determine the region of highest temperature and in order to 
determine the admissible operating temperatures. 

A method for the determination of these temperature distri- 
butions and of the heat-transfer characteristics of a special type 
of heat exchanger is developed. The heat exchanger is composed 
of polygonal flow passages and the passage walls are uniformly 
heated by internal heat sources. The coolant flow within the 
passages is assumed to be turbulent. The circumferential 
variation of the local heat-transfer coefficients is estimated from 
flow measurements made by Nikuradse, postulating similarity 
between velocity and temperature fields. Calculations of tem- 
perature distributions based on these heat-transfer coefficients 
are carried out and results for heat exchangers with triangular 
and rectangular passages are presented. 

INTRODUCTION 

The conventional recuperative type of heat exchanger 
consists of passages for two liquids or gases separated by a 
heating surface. Heat from an outside source is carried with 
a fluid flowing through one of the passages and is transferred 
in the heat exchanger to a second fluid flowing through the 
other passage. Yery often such a heat exchanger is com- 
posed of a large number of tubes, with the two liquids flowing 
inside and over the outside of the tubes, respectively. 

- The regenerative type of heat exchanger has passages for 
one fluid only. During the heating period, heat from an 
outside source is carried to the heat exchanger by a hot fluid 
and is stored within the solid walls of the passages. This 
heat is then given off to a cold fluid, which passes through 
the heat exchanger during the cooling period. 

In this report, a heat exchanger is considered that differs 
from the regenerative type only by the fact that the heat is 
generated by heat sources within the passage walls and is 
transferred to a coolant flowing continuously through the 
passages. The passage walls of this heat exchanger are 
assumed to be flat plates assembled to form a honeycomb; 
thus the flow passages formed in this manner have a polygonal 
cross section. A cross-sectional view of a typical heat 


exchanger of this type is shown in figure 1. The exchanger 
is composed of a number of plates a, which form the coolant 
passages b. The flow of the coolant is normal to the cross 
section shown. High temperatures may be anticipated near 
the comers c of the passages, inasmuch as the rate of heat 
transfer there is expected to be poor. A theoretical inves- 
tigation of the temperature distribution in such a heat 
exchanger was made at the NAOA Lewis laboratory during 
1950 and is presented herein. 

The basis for this investigation is a knowledge of local 
heat-transfer coefficients in passages of noncircular cross 
section. Some information is available on the average 
heat-transfer coefficients in such tubes (references 1 and 2). 
The result of these investigations is essentially that the ex- 
pressions derived for the heat-transfer coefficients in circular 
tubes apply for other cross sections as well, provided the 
diameter is replaced by the hydraulic diameter. (The 
hydraulic diameter is defined as four rimes the cross-sectional 
flow area divided by the circumference of the passage.) 

A knowledge of local heat-transfer coefficients in noncir- 
cular passages is important not only in the present problem, 
but also for several other engineering problems, such as the 
determination of local temperatures in the walls of air-cooled 
turbine blades. No information on local heat-transfer 
coefficients was found in the available literature, however. 
These values are therefore estimated from flow measure- 
ments made by Nikuradse (reference 3) on the basis of the 
similarity between temperature and velocity fields. 



1 Supersedes NACA TN 22S7, “Temperature Distribution. In Internally Heated WbUs of Heat Exchangers Composed of Nonolrcnlar Plow Passages” by E. E. G. Eckert and George M. 
Low, 19J1. 
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SYMBOLS 

WT~ 

The following symbols are used in this report: 

A cross-sectional flow area (sq ft) 

A* ratio of wall area to flow area (dimensionless) 

C internal circumference of passage (ft) 

C* CjD (dimensionless) 

Cf specific heat at constant pressure, Btu/(lb) (°F) 

D hydraulic diameter, 4A/C, (ft) 

g acceleration due to gravity (ft) /(sec 8 ) 

h local heat-transfer coefficient, Btu/(sec) (sq ft) (°F) 

1 ro 

k average heat-transfer coefficient, ^ I h dx, 

Btu/(sec) (sq ft) (°F) 
h* hfh (dimensionless) 

k thermal conductivity of wall material, 

Btu/(sec) (ft) (°F) 

k t thermal conductivity of coolant, Btu/(sec) (ft) (°F) 
k* k/kg (dimensionless) 

A 7 residual value (dimensionless) 

Nu Nusselt number, HD/k c , (dimensionless) 

n coordinate normal to passage wall (ft) 

n* n/D (dimensionless) 

Pr Prandtl number, vjot, (dimensionless) 

p pressure (lb)/(sq ft) 

g local rate of heat transfer, hd, Btu/(sec) (sq ft) 

g* h* d* (dimensionless) 

R radial coordinate (ft) 

Re Reynolds number, uD/v, (dimensionless) 

r rate of internal heat generation, Btu/(sec) (cu ft) 

s wall thickness (ft) 

8* s/D (dimensionless) 

T local total temperature of coolant, °F 

T b bulk total temperature of coolant, °F 

T* Tk/rD 2 (dimensionless) 

t local wall temperature, °F 

t* tk/rD 2 (dimensionless) 

u, v velocity components in x- and y-direction, re- 
spectively (ft/sec) 

x, y, 2 Cartesian coordinates .... J 

x* xjD (dimensionless) 

y* yjD (dimensionless) 

2 * z/D (dimensionless) 

a thermal diffusivity, k g /pgc Pl (sq ft)/(sec) 

/3 angle subtended by two adjacent sides of polygonal 

passage (deg) 

A increment of length (ft) 

A* A/D (dimensionless) 

S increment or difference 

e B turbulent diffusivity of heat (sq ft) /(sec) 

t M turbulent diffusivity of momentum (sq ft) /(sec) 

0 temperature difference, t-T B , °F 

0* dk/rD 2 (dimensionless) 

v kinematic viscosity (sq ft) /(sec) 

p mass density (lb) (sec 2 ) /(ft 4 ) 

local wall shear stress (lb) / (sq ft) 


Subscripts: 

c . conditions for equivalent circular tube 
m .conditions at center of flow passage 

n conditions normal to passage wall 

8 conditions at wall surface 

Superscript: 

* dimensionless quantity 

ASSUMPTIONS 

If this analysis were to be made without any simplifying 
assumptions, the simultaneous solution of the equations of 
motion of the coolant and of the heat flow in the coolant and 
in the passage walls would be required. The following as- 
sumptions are made in order to make these equations 
amenable to solution without seriously curtailing the results 
of the analysis: 

(1) The Prandtl number of the coolant used in the heat 
exchanger is in the neighborhood of 1. This condition is 
well fulfilled by gases and by water above a temperature of 
200° F, excluding the neighborhood of the critical point. 

(2) The passages of the heat exchanger are long enough so 
that in the cross section investigated the flow is fully devel- 
oped, which means that the velocity profile does not change 
its shape in the direction of the tube axis. 

(3) The rate of heat generation in the walls of the heat 
exchanger is uniform. As a consequence of this condition, 
the temperature within the coolant and the walls increases 
linearly in a downstream direction provided the flow is ther- 
mally developed. For a fluid with a Prandtl number of 1, 
the points of thermal and velocity development in a tube 
practically coincide provided that the heating of the tube 
starts at the entrance section. 

(4) The temperature gradient along the tube axis is as- 
sumed small as compared with the gradients in any cross 
section of the passage. 

(5) The thermal conductivity of the solid material is large 
as compared with the thermal conductivity of the coolant. 
This condition is always fulfilled for metal walls regardless 
of the type of coolant used provided that assumption (1) 
applies and for nonmetallic walls if the coolant is a gas. 
For the case of nonmetallic walls and liquid coolants, the 
applicability of the calculations presented in this report must 
be checked in each individual case. Furthermore, the tem- 
perature differences within the passage wall at any one cross 
section are postulated to be small as compared with the 
temperature differences between the wall and the core of 
the coolant. As a consequence of the assumption listed in 
this paragraph, the heat transport within the coolant, normal 
to the tube axis and parallel to the walls, is small as com- 
pared with the heat conduction within the walls and may 
therefore, be neglected. The heat transport within the 
coolant, normal to the tube axis and normal to the walls, 
is of course taken into account. 

(6) The turbulent diffusivity of momentum and that of 
heat t H are equal. 
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FLOW IN TUBES WITH NONCIRCULAR CROSS SECTION 

A thorough, investigation of the flow through tubes with 
noncircular cross sections was made by Nikuradse (refer- 
ence 3). This investigation is used in the present report 
as the basis for estimating local heat-transfer coefficients. 

The flow of water through tubes of several different shapes, 
as shown in figure 2, was investigated in reference 3. Three 
of the tubes had triangular cross sections; an equilateral 
triangle, an isosceles right triangle, and a right triangle with 
the sides enclosing the 90° angle having a ratio of 1 to 2.32. 
One tube had a trapezoidal cross section and two tubes were 
circular with one and two grooves, respectively. In a second 
report, Nikuradse investigated tubes with a rectangular cross 
section; however, only a summary of this report is available 
(reference 4) . The hydraulic diameter of the passages varied 
from 0.3 to 0.6 inch and the length-to-hydraulic-diameter 
ratio varied from 100 to 200. The Reynolds number based 
on the hydraulic diameter and the mean velocity ranged 
from 77,000 to 120,000. Because it is known that the shape 
of a turbulent velocity profile changes only slightly with 
Reynolds number, the results of the calculations should be 
applicable for a fairly large range of Reynolds number in 
the turbulent region. The velocity profiles were measured 






Figure 3.— Velocity contours En triangular passage as measured, in reference 3. Reynolds 

number, 81,000. 


with a small total-head tube in a cross section near the 
downstream end of the tube. 

Results of reference 3 for the isosceles right triangle are 
shown in figures 3 to 5. Similar results for all other cross 



Fkjube 4.— Secondary Sow In triangular passage deduced from measurements of reference 3. 



Fioube 5. — Shear stress distribution on drcumferenee of triangular passage as determined 

in reference 3. 
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sections can be found in the original report. In figure 3 the 
lines of constant velocity are presented as contour lines. 
These contours are indented near the center of each side of 
the passage. The conclusion of reference 3 is that these 
indentations indicate a secondary flow normal to the main 
flow direction. A qualitative sketch of this secondary flow 
is presented in figure 4. Prandtl (reference 5) has ascribed 
the secondary flow to a turbulent mixing motion within the 
fluid, which is more intense in the direction parallel to the 
wall than normal to it. This secondary flow tends to equal- 
ize the velocities and temperatures within any cross section 
of a noncircular tube and is therefore favorable for the simi- 
larity consideration that will be used to deduce the tempera- 
ture field from the measured velocity field. On the other 
hand, because no quantitative knowledge of the secondary 
flow exists, an exact theoretical calculation of the tempera- 
ture field in noncircular passages is impossible. The local 
wall shear stresses were computed by Nikuradse, using the 
measured velocity profiles together with the assumption that 
the Blasius pipe resistance law of the turbulent velocity 
profile for circular tubes applies also for noncircular passages 
on normals to the walls (fig. 5). A check by Nikuradse of 
the calculated average shear stresses against the measured 
pressure drop showed good agreement. 

SIMILARITY BETWEEN VELOCITY AND TEMPERATURE 
PROFILES FOR Pr=l 

It is well known that the velocity and temperature profiles 
in the boundary layer of a fluid with a Prandtl number of 1 
are similar in shape in the absence of a pressure gradient. 
This similarity is immediately apparent from a comparison 
of the momentum and energy equations of the boundary 
layer for two-dimensional steady flow along a flat plate 
(reference 6): 

'(“§ +, 5jD = ^[ ,< ' +, " > l|0 (1> 

p ( u !? + '§D“l;[ p( ’ +, ' ) !tO <2) 

where T represents the total temperature, the Prandtl num- 
ber is equal to 1, and the specific heat is constant. The 
effect of internal friction on the temperature profile is taken 
into account by basing equation (2) on the total temperature. 

Because it is assumed that the turbulent diffusivity of 
momentum e M and of heat e B are equal, equations (1) and 
(2) are similar. The solutions of the equations, namely, 
the velocity and temperature profiles, are therefore also 
similar, provided the boundary conditions are similar. The 
boundary conditions for the velocity field on a flat plate are 
that the velocity is 0 along the wall and has a constant value 
outside the boundary layer. Similar boundary conditions 
for the temperature field are that the temperature is constant 
both along the wall and in the main stream' 

It is therefore evident that similarity between velocity 
and temperature profiles in the boundary layer over a flat 


plate exists in an exact mathematical sense. The same is 
not true, however, for fully developed pipe flow. The bal- 
ance of the forces and of the heat energy on a stationary 
annular volume element with a radius JR, thickness dR, and 
length (in the direction of the tube axis) dz, leads to the 
following equations for fully developed flow in a circular 
tube: 

l^?[ i?p(a+6w) S ==pti |f ^ 

where the heat-conduction term in the 2 -dircction is neglected 
in equation (4). It is apparent that equations (3) and (4) 
are not similar, even if the viscosity and diffusivity terms in 
the two equations are equal. The following type of analysis 
can be made, however: 

The turbulent diffusivity of momentum e. v can be cal- 
culated from equation (3), provided that the velocity profilo 
u=j{R) and the pressure drop dpjdz are experimentally 
known. If it is again assumed that e M and are equal, the 
temperature profile can then be calculated from equation (4). 
Such an analysis, based on constant property values and with 
internal friction neglected, was carried out by Latzko (refer- 
ence 7). The results of this analysis were used to calculate 
the relation between the temperatures and the velocities for 
hydrodynamically and thermally developed flow of a fluid with 
Pr=l. Figure 6 represents a plot of the local temperature- 
difference ratio {T—t,)l{T m —t,) against tho local veloc- 
ity ratio uju m . It can be seen that this relation is very 
nearly linear. This linearity means that the velocity and 
temperature profiles are similar as far as practicality is con- 



Figcbe 8.— EeEMon between local temperatures and local velocities In fully developed tur- 
bulent region of olrcular tube. Based on reference 7. 
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Dimensionless distance from corner, jf-o-/D 

Figure 7 . Local beat-transfer coefficients in comers of non circular passages. 


cemed, eren though, they are not- similar in a strict mathe- 
matical sense. Most of the recent theoretical investigations 
of turbulent heat transfer are therefore based on the assump- 
tion of simil arity between these profiles for a fluid with a 
Prandtl number of 1. 

It is known that the fully developed turbulent velocity 
profile in a circular tube is such that the essential part of the 
velocity change from 0 at the wall to a maximum value at the 
center occurs in a narrow strip around the periphery of the 
tube. The investigations of references 3 and 4 show that this 
condition is also true for tubes of polygonal cross section, as 
long as the included angles between two adjacent passage 
walls are not very small. Nikuradse further shows that the law 
for the velocity variation on a normal to the wall established 
for turbulent flow in circular tubes holds also for the non- 
circular passages in this region of essential velocity variation. 
It can therefore be expected that the similarity between 
temperature and velocity profiles will also be fulfilled reason- 
ably well on normals to the walls for noncircular tubes and 
fluids with a Prandtl number of 1, although a balance of forces 
and energies similar to equations (3) and (4) shows that the 
similarity cannot exist exactly in regions very dose to the 
comers. 

An immediate result of the similarity between velocity and 
temperature profiles is the fact that the local wall shear stress 
determined by the velocity gradient at the wall is propor- 
tional to the local heat-transfer coefficient determined by the 
temperature gradient at the wall. (The same result can be 


obtained from Reynolds’ analogy.) The proportionality of 
these values, together with the knowledge that the average 
heat-transfer coefficients for circular and noncircular tubes 
with the same hydraulic diameter are equal, can be used to 
obtain local heat-transfer coefficients in noncircular tubes 
from the wall-shear-stress data of references 3 and 4. A 
generalization of this result is possible because, for a given 
Reynolds number, the shear stress distribution in the vicinity 
of a comer of the passage depends primarily on the conditions 
near that comer and not on the shape of the passage. 

The local heat-transfer coefficients, as presented in figure 7, 
were obtained by correlating the results of references 3 and 4. 
In this figure the ratio A* of the local heat-transfer coefficient, 
to the mean value is plotted against the dimensionless 
distance from the comer x*. The included angle of two 
adjacent passage walls is the parameter for the curves. Some 
of these corner angles were represented on more than one of 
the cross sections investigated in reference 3, and for these 
angles the curves of A* against x* agree reasonably well. 

The application of figure 7 can best be explained with the 
aid of an example. Suppose it is desired to find the distribu- 
tion of A* over the short leg of an isosceles right triangle, as 
shown in figure 8. The two ends of the curve can imme- 
diately be transposed from figure 7. The remainder of the 
curve is then extr apol ated so that the area under the curve 
divided by length AB is equal to 1. This procedure gives a 
reasonably good approximation as long as the maximum 
value of A* is not considerably greater than 1. 
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DIMENSIONLESS VARIABLES 


With the distribution of the local heat-transfer coefficient 
around the circumference of the passages established, the 
only problem remaining to be solved is the calculation of the 
heat-conduction process ■within the passage walls. Before 
this problem is taken up, an investigation is made to deter- 
mine the dimensionless moduli on which the temperature 
distribution in a heat exchanger of the type under considera- 
tion depends. This investigation may be useful as a basis 
for experimental investigations. None of the simplifying 
assumptions, as summarized previously, is necessary for the 
development in this section of the report. 


The heat conduction within the solid walls of a heat ex- 
changer with internal heat generation is described by Pois- 
son’s equation 


k 



0 


(5) 


The number of parameters for this equation can be reduced 
by the use of the following dimensionless values: 




tk 

rD 2 


X= D 


V *=JL 

y D 


* » 
3 ~D 


( 6 ) 


With these values, equation (5) can be rewritten 


bH* , bH* , m* 
bx* 2± by* 2± bz* 2 


1 = 0 


(7) 


The temperature field that results under the influence of the 
internal heat generation in the solid walls depends on how 
the heat is transferred from the wall surfaces to the coolant. 
The heat must be conducted within the solid material to the 
surface and from there into the coolant. Along the surface, 
therefore, the following boundary condition exists: 



where n represents the direction normal to the surface. In 
terms of dimensionless valuables, this equation becomes 



Figure 8.— Determination of local beat-transfer coefficients In triangular beat-exchanger 

passage. 


The heat transfer within the coolant by conduction and 
convection is governed by an energy equation similar to 
equation (4) and the flow of the coolant is determined by a 
momentum equation (similar to equation (3)) and the cor- 
responding continuity equation. Textbooks on heat transfer 
(for example, reference 8) show that the dimensionless tem- 
perature described by the energy equation of the coolant 
depends, for low-velocity flow and constant property values, 
on the Reynolds number Re and the Prandtl number Pr. At 
high velocities and variable property values, there is an 
additional influence of the Mach number and other dimen- 
sionless expressions characterizing the temperature depend- 
ency on _ the property values. By neglecting the last- 
mentioned influences and summarizing all the factors that 
influence the heat flow in the solid material and in the 
coolant, a functional relation of the following type can be 
deduced : 

d*=f(Re, Pr, k*, x*, y* t z*) (10) 

where the temperature difference 6* is introduced because 
only temperature differentials appear in the equations. The 
function./ as expressed by this equation depends on the 
geometric configuration of the heat-exchanger passages. 

Whenit can be assumed that the lieat-transfor process from 
the walls to the coolant is not influenced by the temperature 
distribution within the wall, the number of factors influenc- 
ing the problem can be considerably reduced. This assump- 
tion is made throughout the calculations in this report and 
seems reasonable as long as previously mentioned assump- 
tion (5) Jholds. In this case equation (8) can bo replaced by 

(11) 

where h is a known function of the flow parameters. By 
introducing the ratio h* of the local heat-transfer coefficient 
to the value averaged over the circumference of the passage 
and the Nusselt number Nu—TiDjk t based on this average 
value, and by changing to the temperature difference 6, 
equation (11) can be transformed to 

(^)r~ Nw ^ 9 * (12) 

The dimensionless temperature within the 6olid walls of the 
heat exchanger is determined by equations (7) and (12) and 
can be presented as a function of the following kind: 

**=/(w«p» **,»*,**) (13) 

This dimensionless temperature difference depends only on 
the dimensionless local coordinates and on the parameter 
h* 

iVup- Temperature distributions and heat-tran9fer char- 
acteristics for any given geometric configuration are therefore 
presented as a one-parameter family of curves. 
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CALCULATION OF TEMPERATUBE DISTRIBUnONS WITHIN 
PASSAGE WALLS 

The problem of calculating temperature distributions within 
the passage walls can be classified in two general categories, 
depending on the relative thickness of these walls. If the 
dimensionless wall thickness s* (see fig. 8) is small as com- 
pared with the length of a wall AB, the temperature differ- 
ences in the direction normal to the wall surface are small as 
compared with the temperature difference in the direction x* 
parallel to the surface. Only the temperature differences in 
this parallel direction need therefore be considered. Here- 
inafter this special case of the problem is referred to as the 
“one-dimensional problem.” On the other hand, if s* be- 
comes large, the heat flow in both the z * and s * direc tions must 
be considered. This second and more general problem is 
referred to as the “two-dimensional problem.” 

One-dimensional problem. — Consider a thin plate of 
thickness s that is a section of the heat-exchanger passage 
walls. Heat is generated uniformly throughout the plate 
at a rate r. The local plate (or wall) temperature t is assumed 
to be varying in the z-direction only. The bulk temperature 
of the coolant is T s and the thermal conductivities of the wall 
and coolant are k and k e , respectively. The local coefficients 
of heat transfer to the coolant above and below the plate are 
denoted by h and h', as shown in the following sketch: 


h n Q a 



d*0* A*+A”\ r ^ , n 
^ i -^ F ^Nu8*+ 1 =0 (16) 

This equation expresses the dimensionless wall temperature 
6* as a function of the dim ensionless coordinate x * and a 
single parameter ( h*+h'*)Nufs*k *. For any given geo- 
metric configuration (which determines A* and k'* as func- 
tions of x*), the dimensionless wall temperature is a function 
only of the average Nusselt number Nu, the dimensionless 
wall thickness s*, and the ratio of the thermal conduc- 
tivities k*. 

In general, equation (16) must be solved by numerical 
means inasmuch as A* and h'* are experimentally determined 
functions of x*. Several numerical methods of solution can 
be applied and two of these are considered here. The relaxa- 
tion method (see, for example, reference 2, pp. 365-379, or 
reference 9) has the advantage that it is easy to apply and 
that computational errors are immediately apparent. A 
second method of solution, based on the Runge-Kutta 
method (reference 10) is presented in the appendix. This 
method is preferable when results of high accuracy are 
desired. It is less convenient, however, than the relaxation 
method because computational errors are much more difficult 
to detect. 

For the purpose of solving equation (16) by the relaxation 
method, the equation is first expressed in finite-difference 
form. Consider a grid, or net of points, placed into the wall, 
any two adjacent points being separated by a small but 
finite distance A* as indicated in the following sketch: 


m 



n 



Hi 



m 


■■■■ 




The heat balance for an element of volume with the 
dimensions dx and a and of unit depth is 

QA-Q^Q' A-Qa^rQa (14) 

where 

e , =- i »(s+S^) 

Q t =rsdx 

Qa-\- Qa=(h + h r ) (f — T b ) dx 

"With the preceding values of Q, equation (14) becomes 


The second derivative of 0* at an arbitrary point 0 can be 
expressed in terms of the temperatures 8* at this point and 
the two adjacent points as follows: 

/d$*\ (dB*\ .. ft*— 

i\dz*/ 0 ~ A* a\dx*) 2 ~ A* 

d* e* i r /d0*\ (dQ * M 

dx* 2 ~A* Li\dx*J 0 o\dx*JJ 

«^(C+e 2 *-2e 0 *) (17) 

With this value for the second • derivative, equation (16) 
becomes 


dH h+h' 
dx 2 ks 


(f-r s )+£=o 


(15) 


0L*+ e t *~ 00* [2 +^(A*+A'*)oA«]+A«=0 (18) 


Equation (15) can be expressed in terms of dimensionless 
variables as 


The relaxation method, used for solving this equation, is 
subsequently discussed. 
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Two-dimensional problem. — Consider a slab of homogene- 
ous solid material with unit thickness, as shown in the 
following diagram : 



Heat is again generated within the material at a uniform 
rate r. The heat balance for an element of volume 'with the 
dimensions dx and dy and of unit depth is 

&+&+&=&'+&' (19) 

A constant flow of heat through the slab normal to the 
plane of the preceding sketch does hot influence this heat 
balance. Local deviations from this constant flow of heat 
are neglected. 

The individual terms in equation (19) Are 


where n is the direction normal to the surface.. Equation (22) 
can again be written in terms of the dimensionless variables 

Equation" (21) together with boundary condition (23) fully 
describes the two-dimensional problem. All physical vari- 
ables are.again grouped into a single parameter h* Nu/k*, 
which appears in the boundary condition. In order to apply 
the relaxation method of solution to these expressions, they 
are first converted to difference equations. 

For this purpose a rectangular grid is placed into the slat), 
as indicated in the following diagram: 



Q,=—kdx-^ 

q,‘--kix(~+~^Av) 

Qi=rdxdy 

Equation (19) can' therefore be written 

bH . bH _ r 
bx 2 'by 2 ~ k 

Or, in terms of the dimensionless variables, 

b 2 9* &e* 
bx* 2+ by* 2 1 

If the slab is bounded by a coolant whose bulk temperature 
is T b and if the local surface heat-transfer coefficient at a 
point s on the surface is h, the boundary condition at that 
point is .. 

Kl).— 


Adjacent net points are separated by a distance A*. Tho 
derivatives of 9* at an arbitrary point 0 can be expressed in 
terms of the temperature function at surrounding points. 
Thus the first derivatives are 



a * a * 



bx*J 

1 ~ A* 

0 “ 

o\bx*J 

A* 

b9*\ 

. . 9*— 9* 

(K) 

a * a ♦ 

1 « 0 


'o A* 

0 W) 

« A* 


and the second derivatives become 

b*d* 0i*+0a* — 20 o * 

5^2“ A* 2 

b 2 9* 6 t *+9 i *-29 0 * 

dy* 2 ” A* 2 

With these values, equation (21) becomes 

0i * + 0 2 * + 0,* + 0<* - 4 0„* + A**= 0 (24) 

The boundary condition can be evaluated by assuming that 
points 1, 0, and 3 lie on the surface of tho wall. Point 4 then 
lies in the stream and its temperature must bo expressed by 
the normal derivative at point 0. 

/50*\ 9* — 9* 

W/ o 


( 20 ) 


( 21 ) 


( 22 ) 
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or, with the use of the value of <56* [dy* as given by equation 
(23), 

0 4 *=0q*( (25) 

Along the boundary, therefore, the following equation 
applies: 

eS+es+ef- e 0 * (s +^i- a*)+a* 2 =o (2 6) 

Equation (24) together with boundary condition (26) can 
again be solved by means of the relaxation method. 

Solution by relaxation method. — A heat exchanger com- 
posed of a large number of rectangular passages will be dis- 
cussed. The walls of the heat exchanger are made of a 
homogeneous material in which heat is generated at a uni- 
form rate r. Figure 9 represents a cut through the heat 
exchanger so that the flow of coolant is in a direction normal 
to the cut. In the discussion that follows, it is assumed that 
the geometry of the configuration is given in the dimensionless 
system of coordinates. 

A complete discussion of the relaxation method is not 
presented herein, inasmuch as it is generally available 
elsewhere. (See, for example, reference 9.) The essential 
features of the method can be outlined as follows: 

Suppose it is desired to solve a given finite-difference 
equation over a certain area of integration. The equation 
is of the following type: 

/(fl*)=0 (29) 

The solution of equation (29) must also satisfy prescribed 
conditions at the boundary of the area of integration. 

First, it is necessary to select a number of net points cover- 
ing the entire area of integration. The distance A* between 
net points is arbitrary, with the accuracy of the final solution 
increasing as the distance between points is decreased. 
Next, values of 6* are assumed at each net point. If by 
chance these assumed values of 6* are the correct values, 
then they satisfy the appropriate finite-difference equation 
at all net points. 

In general, however, the assumed values of the function do 
not satisfy the difference equation and the left side of equa- 
tion (29) is equal to some residual value N instead of zero. 
At any given net point 6* must then be adjusted in order to 
make N vanish at that point. This adjustment of 9* also 
changes the residuals at adjacent net points. However, 
if this process of adjustment is started at the point at which 
the absolute value of N is greatest and is then repeated for 
points at which the value of the residual is successively less, 
the correct values of 9* for the entire net eventually are 
obtained. 

In applying this method to the heat exchanger under 
consideration, it is first assumed that the wall thickness s* 
is small, so that the one-dimensional solution applies. 



FinrEE 10. — TUn-waHed heat exchanger with rectangular passages. (Calculation Is confined 

to region ABOJ 



FictDXE U.— Relaxation net lor rectangular thin-walled heat exchanger. 


The actual configuration to be discussed is shown in fig- 
ure 10. Because points A, B, and C are points of sym- 
metry, only the region bounded by these points need be 
considered. Furthermore, the distribution of temperature is 
also symmetric about point O with respect to AO and 00. 
An enlarged view of the section under consideration is 
represented in figure 11. 

The section has been subdivided into a number of net 
points a, b, c , . . . h, each separated by a distance A* from 
adjacent points. The first step in the analysis is the assump- 
tion of temperatures at all net points, which can be done by 
setting up a balance between the total heat produced within 
the walls and the total outflow of heat from the walls. In 
terms of the dimensionless variables, this heat balance 
becomes 

(j* pj £* dx*= 1 (30) 


If, for the ini ti al assumption, h*=h'*—l and 6* is constant, 
the following expression is obtained: 



(31) 


The value of 9* given by equation (31) is assumed to exist 
at all net points. The appropriate difference equation at 
points 6, c, d,f, and g is (from equation (18)) 
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0i* + 0 a * - 0o* [2 ■ + (A* + A'*)«A* 3 ] + A «= No (32) 

where the subscript 0 refers to - the point at which the equa- 
tion is to be applied and subscripts 1 and 2 refer to adjacent 
points. At point a, which is influenced by three points, the 
following equation is required: 

2 6* + 9* — 0„* [ 3 + -pjx (A* + A/*)oA* J ] + A* 3 = iVa (33) 

At point e, which is a point of symmetry, the following 
expression applies: 

20/-0.*[2+^(A*+A'*) 4 A* , ]+A* s =fV < (34) 

A similar expression applies at point A. With these equa- 
tions, the residuals N are calculated at each net point. The 
point at which N has the largest absolute value is then 
selected and 0* at that point is adjusted so that the residual 
vanishes. The effect of this adjustment on the residuals at 
adjacent points is calculated with the aid of the appropriate 
finite-difference equation. 

The process is repeated for the point at which the next 
largest value of the residual appears. Eventually, if the. 
process is repeated often enough, the residuals at all net 
points approach zero. The final adjusted values of 6* then 
satisfy the appropriate equations at all net points. 

These values of 0* should also satisfy the heat balance as 
given by equation (30). This equation can therefore be 
used to check the validity of the final temperature distribu- 
tion. If the check is unsatisfactory, it is necessary to 
continue the solution of the finite-difference equations by 
using a smaller net spacing. 

The same heat-exchanger configuration is now investigated 
without making the assumption that the wall thickness is 
small. The two-dimensional equations are therefore applied. 
An enlarged view of the shaded portion of figure 9 is shown 
in figure 12. The lines CD, FG, and GA are lines of sym- 
metry, and there is no flow of heat acros s these lines. The 
temperature distribution along line DEF will be symmetric 
about point E, but there will be a flow of heat across this 
line. A large number of net points are selected to cover the 
entire section of the configuration shown in figure 12. (Only 



£ D 

Figues 13.— Relaxation net tor rectangular thlok-walled beat exchanger. 


a few ol these points are presented in the fig.) An initial 
value for the temperature along lines AB and BC is obtained 
by setting up a balance between the total heat produced 
within the section and the flow of heat to the gas. In 
terms of the dimensionless variables, this heat balance 
becomes 

<kf. (35 > 

If, for the initial assumption, A*=l and 0* is constant, 
then the following expression for 0* is obtained: 


° ~ 4 Nu 


It can now be assumed that the temperature as given by 
equation (36) exists along the surface of the wall and that 
a somewhat higher temperature exists at internal points. 
The assumed temperatures are then adjusted by means of 
the relaxation method. 

The appropriate finite-difference equation for internal 
points, such as point e, is written as follows: 

06* + 0* 4- 0** + 0 tf * — 4 0, + A*’ = JV, (37) 


For points along the boundary, such as point b, the equation 
is 



and for points on a line of symmetry, such as point/, 

0 c *+20/+0 t *-40/+A **=1V) (39) 

The appropriate equation for points along DEF can be 
determined by taking advantage, of the antisymmetry about 
this line. At point k, for instance, the equation becomes 

0**+0»*+0n*+0/-40 t *+A* 2 =jV t (40) 


Expressions of the type (37) to (40) apply at all net points. 
The method of determining the actual temperature distribu- 
tion is the same as the method outlined in the previous 
section. 

The final temperatures along tho surface of the wall can 
be checked with the aid of equation (35). If the check 
is unsatisfactory, a finer net spacing is required. 

RESULTS 

One-dimensional solution. — The ono-dimensional solution 
was applied to heat exchangers composed of rectangular and 
triangular passages, respectively. In each case the passages 
were staggered in order to minimize tho expected hot spots. 

The rectangular configuration is represented by figure 10. 
The height-to-width ratio of each passage in this configura- 
tion is 1 to 5. Figure 1 represents the triangular configura- 
tion. Each passage in this configuration is an isosceles 
right triangle. 
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Length ratio, x'-x/D 

Fiqcee 13.— Temperature dfetrfbatton for passage of thto-walled triangular beat exchanger 

(flg- 1 ). 

It has previously been stated that, for each configuration, 
the temperature distribution in the passage walls is a function 
only of the parameter Nu{s*k*. Accordingly, temperature 
distributions were calculated for several values of this param- 
eter, ranging from 0.1 to 100. This range of values is 
believed to include all values actually encountered. ■ 

Temperature distributions were calculated according to the 
method outlined in the previous section of this report and 
were then checked by equation (30). Local rates of heat 
transfer g *=fl*0* were calculated for this purpose. Accord- 
ing to equation (30), the mean value of g* when multiplied 
by 2Nu/k*s* should equal 1. Because the finite-difference 
method of solution is essentially an approximate method, the 
results were not expected to be exact. All results presented 
in this section of the report, however, were held to an error 
of less than 6 percent. The final temperatures were multi- 
plied by a constant scale factor -in order to satisfy equation 
(30) exactly. 

Temperature distributions and local rates of heat transfer 
are presented for the triangular heat exchanger in figures 13 
and 14 and for the rectangular heat exchanger in figures 15 


and 16. In the temperature plots (figs. 13 and 15), the 
temperature difference (t — T B ) is referred to the temperature 
differ enc e (t—T B ) c for a thin- walled circular tube with the 
same hydraulic diameter and an internally heated wall with 
the same physical properties. The temperature difference 
for the circular tube is given in the nondimensionat form by 
equation (31). Similarly, the ordinate in figures 14 and 16 
refers the rate of heat transfer g to the rate of heat transfer 
from the wall of a thin circular tube g„. It is believed that 
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Fiouhe I*.— Temperature distribution for passage of thin-walled rectangular heat exchanger 

(flg. 10). 


the interpretation of the temperature and heat-transfer 
curves is simplified by this final change of ordinates. Tem- 
perature distributions for the triangular configuration are 
shown in figure 13. These curves are presented for values 
of Nu/s*k* of 0.1, 1.0, 10, and 100. As was expected, the 
temperatures peak near the comers of the passage. The 
largest temperature differences are encountered for large 
values of the parameter Nu/8*k*. A large value of this 
parameter indicates a low conductivity of the wall material 
or a small wall thickness. Both of these conditions are 
conducive to low rates of heat exchange within the wall. 
Local rates of heat transfer corresponding to the temperature 
distributions just discussed are presented in figure 14. 

Temperature and heat-transfer curves for the rectangular 
configuration are shown in figures 15 and 16. The results for 
this configuration are similar to those for the triangular 
configuration. 

In order to determine when the one-dimensional solution 
can be used, it is necessary to obtain some information on 
the temperature difference that exists across the passage 
walls on normals to the surfaces and to compare this tem- 
perature difference with the temperature differences along 
the surfaces. : An estimate of the temperature difference 
8t n across the wall may be obtained by calculating this value 
for a flat plate. The result of this calculation is 

82 < 41 > 

The equation can be changed to the dimensionless values 
to yield 

5t ** = T " ' ' ' < 42 > 


&t a s*A T u 
0 t-T B ) e ~ 4 k* 


( 43 ) 


The one-dimensional solution applies as long as this value is 
small as compared with the temperature differences presented 
in figures 13 and 15. 

Two-dimensional solution. — The two-dimensional solution 
for the temperature field within the heat exchanger depends 
on two parameters; namely, the dimensionless wall thickness 
s* and the value Nu h*lk*. In addition, the time required 
to obtain the solution for a special case is much longer than 
for the one-dimensional solution. Only one example was 
therefore calculated; namely, the temperature distribution 
within the walls of a heat exchanger composed of rectangular 
passages, as shown in figure 9. The ratio of the two sido 
lengths.of the rectangle is 1 to 5. The ratio of the hydraulic 
diameter to the short side of this passage is 1.67 and the 
dimensionless wall thickness is 0.6. Figure 17 presents the 
results of the calculation using the relaxation method with a 
network of 88 points. Lines of constant temperature (iso- 
therms) are shown within the portion of the heat-exchanger 
walls that is shaded in figure 9. As may be seen, the heat 
flow within this wall is mainly in the direction normal to 
the wall surfaces. The maximum temperature differences on 
any normal to the surface are not very different from the 
value in a flat plate as calculated in equation (42). The 
temperature differences along the surface of the wall are 



This temperature difference can again be referred to the 
temperature difference ( t—T B ) c as follows: 


Figube 16. — Local rate of heat transfer from passage of thln-walled rectangular heat exchanger 

(flg. 10). 
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Fioure 17. — -Isotherms la internally heated walls of thlek-walled rectangnlar heat exchanger. 
(Shaded portion. of fig. 9 shown.) 2Vo/i*fc*, 10.0; ratio of fiow area to toal area, (UlfT, 
»\ 0.8 ;C/i',1|5. 


appreciably smaller than the temperature differences across 
the wall. The values of the one-dimensional solution cannot 
therefore be expected to apply to this case. Actually, the 
temperature difference found in figure 17 along the surface 
of the wall is considerably greater than for the corresponding 
one-dimensional case. In addition, heat is_ generated within 
the comer area of the rectangular-wall configuration and has 
to be conducted away along comparatively long paths. An 
additional temperature increase can. therefore be found 
within this comer area. 

CONCLUSIONS 

A method to calculate the temperature distribution in a 
heat exchanger composed of noncircular flow passages with 
internally heated walls has been presented. 

Local heat-transfer coefficients along the circu mf erence of 
the heat-exchanger passages were obtained from flow 
measurements made by Nikuradse, assuming similarity be- 
tween the velocity and temperature fields. The heat-transfer 
coefficients, as determined in this maimer, decrease sharply 
near the comers of the passages and vanish at the comem. 
This decrease becomes more pronounced as the comer angle 


is increased. Near the comers these coefficients are es- 
sentially influenced only by the magnitude of the included 
angle of the comer. 

It was shown that the dimensionless temperature distribu- 
tion within the passage walls depends on a single parameter, 
provided the dimensionless wall thickness is small. Numeri- 
cal evaluations for triangular- and rectangular-passage con- 
figurations for a wide range of the aforementioned parameter 
were carried out. The results of these evaluations are pre- 
sented so that temperature differences arising in the walls of 
heat exchangers with the investigated passage shapes for any 
condition within the range of practical interest can be read 
off one of several curves. These curves show a temperature 
increase near the comers of the passages. This increase 
becomes more pronounced for high Nusselt numbers based 
on the average heat-transfer coefficient, for low wall-thickness 
ratios, and for low ratios of the conductivity of the wall to 
the conductivity of the coolant. 

The dimensionless temperature distribution for thick- 
waned heat exchangers depends on the wall-thickness ratio 
in addition to the property parameter. Because the time 
required to calculate temperature distributions in the whole 
field of interest determined by the two parameters is prohibi- 
tive, only a specific example was evaluated numerically. 
The method of calculation is presented in great detail, how- 
ever, so that evaluations for other interesting cases can be 
carried out. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, October 4 , 1950 


APPENDIX 

SOLUTION BY RUNGE-KUTTA METHOD 


Equation (16) is of the following form: 

g* +1=0 (Al) 

The boundary conditions are 



where a and b are points of symmetry. 

In general, a solution of equation (Al) could he obtained 
by assuming an initial value of 0* at point a In addition to 
the first of the boundary conditions (A2) and by working the 
solution toward point b using a numerical method of inte- 
gration. If the boundary condition at point b is not satisfied, 


the solution has to be repeated with a new initial value. 
This trial-and-error process is tedious and can be avoided 
by splitting equation (Al) in such a manner as to keep an 
undetermined constant c in the solution. This constant is 
finally determined by satisfying the second boundary con- 
dition. Let 

e*=e l *+cd 2 * (A3) 

0i*+ 1=0 

2S-/CO 02 *= 0 

where 


e *( a )= 0 0 3 *(a)=l 
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These expressions still satisfy equation (Al) and the first 
of the boundary conditions (A2). After the solutions for 
and 0 2 * have been numerically obtained, the constant c 
is determined so that the boundary condition at point b is 
fulfilled; thus 



The final temperature distribution is given by equation (A3). 

The temperature distribution in a thin-walled heat ex- 
changer with triangular passages was calculated in this 
manner (for A%/s*&*=10) and was compared with the cor- 
responding curve in figure 13. The Runge-Kutfca method 
(reference 10) was used for the numerical integration. This 
particular comparison showed differences of 2.2 percent of 
the temperature ratio near the comers of the passages and 
smaller deviations elsewhere. The temperature curves are 
uncertain to the same order of magnitude, however, because 


of the freedom in the extrapolation of the heat-transfer- 
eoefficient curves. The simpler relaxation method is there- 
fore regarded satisfactory for the present purpose. 
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